Octahedral Engineering of Orbital Polarizations in Charge Transfer Oxides 
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Negative charge transfer ABO3 oxides may undergo electronic metal-insulator transitions (MIT) 
concomitant with a dilation and contraction of nearly rigid octahedra. On both sides of the MIT are 
in-phase or out-of-phase (or both) rotations of adjacent octahedra that buckle the B-O-B bond 
angle away from 180°. Using density functional theory with the PBEsol+C/ approach, we describe 
a novel octahedral engineering avenue to control the B 3d and O 2p orbital polarization through 
enhancement of the BOe rotation "sense" rather than solely through conventional changes to the 
B-O bond lengths, i.e. crystal field distortions. Using CaFeOa as a prototypical material, we show 
the flavor of the octahedral rotation pattern when combined with strain-rotation coupling and thin 
film engineering strategies offers a promising avenue to fine tune orbital polarizations near electronic 
phase boundaries. 

PACS numbers: 75.25.Dk, 61.50.-f. 71.20.-b 



I. INTRODUCION 

Transition metal oxide (TMO) perovskites are known 
to be strongly correlated materials,^ whose properties 
are controlled by a complex interplay between geometric 
and electronic degrees of freedom. These are determined 
by considering the relative magnitude between various 
energy scales and interactions: the energy difference of 
the transition metal (M) d orbitals and the oxygen p 
states, referred to as the charge transfer energy, and the 
strength of the on-site Hubbard U interaction of the d 
electrons. The charge transfer energy is more important 
in TMO where low-energy excitations are of p— d-type, 
whereas the Coulombic interaction, which localizes the 
electrons on the M-site, produces the insulating state 
in Mott-Hubbard systems." The properties of correlated 
electrons of TMO are controlled in part by the relative oc- 
cupancy of the different transition-metal d orbitals.' The 
relative d orbital occupancy is largely determined by the 
crystal field experienced by the transition-metal cation; 
this electrostatic field is the result of the 2p electronic 
density of the coordinating oxygen ligands. The latter, in 
turn, is directed by the extended geometric arrangement 
of the nearest neighboring oxygen atoms. 

Orbital occupancy can be tuned, for example, through 
chemical substitution,"^'' epitaxial strain,''"' or by super- 
lattice formation in thin films. *^^" Isovalent substitutions 
are important to charge transfer-type oxides, because 
rather than modifying the M-site electronic configura- 
tion, cations with different ionic sizes but the same formal 
valence renormalize the transfer interaction and the one- 
electron bandwidths through changes in the crystal struc- 
ture, i.e., interatomic bond angles and distances. The 
crucial distortion in ABO^ perovskites is the buckling of 
the B-O-B bond angles, because the effective d electron 
transfer interaction between the neighboring transition 
metal sites is mediated by the angular overlap with the O 
2p states. ^^ When the B-O-B bond angle deviates from 
the ideal value of 180°, the transfer interaction weak- 
ens and the bandwidth narrows. Such distortions to the 



inter-octahedral bond angles are typical in GdFeOa-type 
perovskites, and therefore knowing how to control the 
amplitude of the bond angles distortions is critical to 
tailoring the charge transfer and d-orbital polarization. 
While it is well-established that greater octahedral rota- 
tions produce more buckled bond angles and narrower 
bandwidth, '^^ it is not well- understood how the rotation 
"sense" — be it in-phase or out-of-phase along a specific 
Cartesian direction — influences the differential charge 
occupancy on the _B-sitc d and O 2p orbitals. 

In this work, we use density functional calculations to 
show that controlling the subtle flavor of the octahedral 
rotation sense is as important as the overall amplitude of 
the rotations when engineering the electronic structure, 
vis-d-vis orbital polarization, of charge transfer oxides 
near electronic MIT. We explain this behavior using the 
definition of orbital polarization V of mn orbital relative 
to 771(2 orbital, 
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where ni-^mn and ni^m^^ are the occupancies of Ihrnn) 
and |^2'77;2) orbitals, with orbital quantum number li 
and magnetic quantum number mn, respectively.'^ With 
this definition, orbital polarization becomes an effective 
measure of the charge excess in the former orbital with 
respect to the latter. The metric then enables us to isolate 
the contribution of each rotation sense to the orbital 
polarization by judicious choice of structural distortions 
and subsequent calculation of V. 

We demonstrate this utility for the specific case of 
the Fe 3d (cg-symmetry) and O 2p orbitals in CaFeOa. 
Specifically, we show that out-of-phase octahedral rota- 
tions drive a transfer of electronic charge to the apical 2p 
orbitals along the Fe-0 chains coUinear with the axes of 
rotation, whereas the Fe Cg state are less sensitive. We 
then use tensile epitaxial strain as a practical handle to 
enhance the orbital polarization using the P-contribution 
from the out-of-phase octahedral rotations to tune the 
ligand charge density — a crucial parameter involved in 



oxide-based MIT materials. Our results indicate that 
electronic-structure engineering requires both the ampli- 
tude and sense of the octahedral rotations be considered 
on equal footing when designing perovskite TMO for 
integration into next-generation electronics.^' 



II. STRUCTURAL AND COMPUTATIONAL 
DETAILS 

To explore the rotation-orbital polarization interac- 
tions, we make use of the prototypical charge trans- 
fer oxide CaFeOg (CFO), because it lacks Jahn- Teller 
distortions^'''^'' and exhibits a sequence of structural 
phase transitions concomitant with a first-order MIT 
near room temperature (RT).'^'''^' Upon cooling below 
290 K, the metallic paramagnet becomes semiconducting 
(space group P2i/n)/''' whereby the appearance of an 
octahedral breathing distortion (BD) together with the 
octahedral rotations, acts to open the electronic gap,''' 
as it has also been observed in rare-earth nickelates.^" 

Besides the octahedral breathing distortion, which con- 
sists of uniform elongation and contraction of Fe-0 bond 
lengths of adjacent octahedra, CaFeOa adopts-^^ the most 
common a^a'^c^ Glazer tilt pattern, ■^^ with buckled Fe- 
O-Fe bond angles (Fig. 1). Glazer notation allows us 
to distinguish the relative sense of the octahedral rota- 
tions about Cartesian axes passing through the Fe cation. 
The -I- (— ) superscript indicates in-phase (out-of-phase) 
rotations about each x-, y- and z-axis. This rotation 
pattern consists of in-phase rotations about the z-axis 
(a'^a'^c^) and is accompanied by two out-of-phase rota- 
tions, a~a'^c^ and a'^a^c'^, of the same magnitude about 
the X- and y-axes, respectively. 

For a perovskite with a Bravais lattice that has orthog- 
onal interaxial angles and the a^a^c^ pattern, the two 
important Fe-O-Fe bond angles are 6xy and 0z (Fig. 1). 
When only in-phase a^a^c'^ or out-of-phase a^a°c~ rota- 



tion patterns are present, the O^y Fe-O-Fe angle measured 
within two subsequent Fe-0 layers orthogonal to the z- 
direction is the same. The phase, i.e., the "rotation sense," 
does not change the rotation magnitude across the layers. 
Naively, one would deduce that the electronic structure, 
and in turn the orbital polarizations, is insensitive to the 
sense of the rotations, since both rotations induce the 
same rotation magnitude about the z-axis — we investigate 
this supposition explicitly. 

To evaluate how the rotation sense alters the electronic 
distribution, we carry out density functional calculations 
using the projector-augmented wave (PAW) formalism''^ 
as implemented in the Vienna Ab initio Simulation Pack- 
age (VASP)^'' within the revised Perdew-Burke-Ernzerhof 
(PBE) generalized gradient approximation for densely 
packed solids'' plus Hubbard U method^'' (PBEsol-|-C/) 
based on the electronic and atomic structural analysis 
reported in our previous study. ^" We chose the spheri- 
cally averaged form of the rotationally invariant effective 
U parameter of Dudarev et al.-^ with a Ucff — 3.0 eV, 
hereafter U, on the Fe d orbitals. Note that no significant 
differences in the main features of the MIT occur for val- 
ues up to U values of 4.0 eV, as previously reported in 
literature, '^^ and we treat the double- counting term within 
the fully localized limit. ■^^ We impose FM order on all 
Fe-sites then fully relax the spin density. We use a plane 
wave energy cutoff of 600 eV, which is then increased 
up to 850 eV for calculations of the electronic densities 
of states (DOS). Ground state structures (DOS) are ob- 
tained by sampling the Brillouin zone with a minimum of 
a7x7x7 (9x9x7) /c-point mesh and the integrations 
performed with 20 meV Gaussian smearing (tetrahedron 
method). Atomic position relaxation are stopped when 
Hellmann-Feynman forces are minimized to a 0.5 meV 
A^^ tolerance. 



III. RESULTS AND DISCUSSION 
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FIG. 1. Rotationally distorted CaFeOg structure with the 8d 
(equatorial) and 4c (apical) oxygen Wyckofl positions indicated 
by blue and red spheres, respectively. Pseudo-cubic axes (a;, y, 
and z) and Fe-O-Fe bond angles Oxy and 6^ are also shown. 



Our starting point is the analysis of the electronic 
DOS of cubic CaFeOg without any octahedral rotations 
[Fig. 2(a)]. The features of interest are in the energy range 
-8 to 3 eV. The O 2p states dominate the entire range, 
validating the charge transfer description. An energy gap 
appears at the Fermi level (Ep) in the spin-down channel. 
The Fe 3d orbitals contribute to the spin-down states with 
a main band centered near -5 eV. In the spin-up channel, 
the Fe 3d states are divided into two principal bands 
centered about -6.6 and -2.2 eV, respectively, contributing 
to the metallicity with a broad band extending up to 2 
eV above Ep. 

We now examine the symmetry reduction consequences 
from the octahedral rotations and breathing distortions 
(BD) on the electronic structure of the cubic phase. The 
BD mode is responsible for splitting of the la Fe site in the 
undistorted structure (Pm3m) into two Wyckoff positions, 
26 and 2c, each centered in smaller and larger octahedra 
arranged as in a 3D checkerboard (Fig. 1). Both the in- 
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FIG. 2. Total and atom-resolved CaFeOa density of states 
(DOS) for structures with (a) Pm3m, (b) PA/nmc and (c) 
lA/m symmetry. The structures in (b) and (c) are obtained 
by applying a breathing distortion together with aPa^c^ and 
aPaPc~ rotations to the cubic structure, respectively, at identi- 
cal breathing amplitudes and d^y rotation angles but different 
rotation sense (cf. filled circles in Fig. 4). Differences at the 
Fermi level (0.0 eV) induced by the two rotation senses are 
nearly indistinguishable. 



phase aPa^c^ and the out-of phase aPaPc" rotations split 
the 3d oxygen site symmetries in the parent phase into two 
distinct sets; the consequence of these distortions is that 
the Fe d-orbitals experience an electrostatic interaction 
from the oxygen ligands which no longer has the complete 
octahedral cubic symmetry. Coupling the BD with in- 
phase or out-of-phase rotations, results in structures with 
PA/nmc or lA/m tetragonal symmetry, respectively. This 
symmetry lowering changes the DOS with respect to the 
cubic phase mainly in energy windows far from Ep: The 
spin-down bandwidth decreases, enhancing the band gap, 
and the spin-up states shift to higher energies in the 
conduction band (near 1 eV). Comparison of the DOS for 
the different rotation flavors [Fig. 2(b) and (c)], however, 
shows no appreciable differences near Ep; the low energy 
bands appear insensitive to the rotation sense. 



A. Orbital Polarization 

We now track the change in the orbital polarization, V, 
as we increase the breathing distortion together with the 
rotation amplitude of the oxygen displacement patterns 
associated with either the a'^a'^c'^ or a'^a'^c" rotations 
[Fig. 3(a) and Fig. 3(b)] about the z-direction. To com- 
pute V, occupancies are calculated by integrating the 
angular momentum resolved DOS for particular \lmi), 
obtained by projecting the Kohn-Sham wave functions 
onto the spherical harmonics, starting from the lower edge 
of the anti-bonding valence band states up to Ep. 

The orbital polarization for the Fe Cg states (/ = 2), 
e.g., Vd2,do which specifies the filling of the da,2_y2 orbital 
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FIG. 3. (a) a^a'^c^ and (b) a°a''c~ rotation patterns of two 
adjacent octahedra along the z-axis. (c) The out-of-phase 
rotations induce a transfer of charge from the equatorial plane 
(shaded), containing the Px and Py orbitals of the apical oxygen 
atom, to the Fe-0 apical bond about the z-axis, increasing the 
electron occupancy of the pz orbital lying along it. (d) On the 
Fe site, irrespective of the rotation sense, octahedral rotations 
induce charge transfer from the d^2_y2 to the d3z2_r2 orbital, 
the latter directed along the apical Fe-0 bond. 



orbital (to; — 0) is 
j;._y2 + 77-3^2 _r2), where 
_j,2 and d^z^-r^ orbital 



(mi — 2) relative to the d^^-. 
'Pd2,do = (rix^^y^ - n^z-^^r^)/{n 
nx2_y2 and n^z'^-r^ are the d; 
occupancies, respectively. We first analyze the orbital 
polarization of the apical [Fig. 4(a-d)] and equatorial 
oxygen atoms in CaFe03 (see Fig. 1) for all possible 2p 
orbital pairs with increasing amplitude of rotation and 
BD under a fixed volume constraint: 'Ppx,py, 'Ppx,pz and 
Vpy^pz {nil = -1,0,1). Next, we evaluate 7'd2,do[Fe(l)] 
and 7^(22,^0 [Fe(2)], i.e., the orbital polarization of each 
unique Fe site. 

The Vpx,py polarization of the apical oxygen atoms 
evolves identically for both in-phase and out-of-phase rota- 
tions (data not shown); the small polarization value (0.1%) 
also indicates that both p^ and Py orbitals are equally pop- 
ulated, consistent with the tetragonal constraints imposed 
by the Bravais lattice. Over the structure-distortion space 
examined,^'' we find the Vpx,pz and Vpy^pz polarization 
is positive with px and py preferentially occupied over 
the orbital directed along z {pz). Both polarizations are 
~ 10% and the width of variation is 4%. These axial 
polarizations diminish for both rotation flavors if the BD 
increases. We understand this behavior as follows: The 
Pz orbital pointing along Fe-0 bond axis increases its 
electron occupancy owing to a charge transfer from the 
equatorial plane, containing px and Py orbitals, towards 
the Fe-O-Fe bond axis [Fig. 3(c)]. 

The orbital polarization involving pz is nearly constant 
with respect to increasing in-phase rotation and turns out 
to have a minimum at large out-of-phase angles and large 
breathing distortion values. At fixed BD amplitude, large 
out-of-phase angles direct the charge transfer towards the 
Fe-0 bond along the (rotation) z-axis. 

For the equatorial oxygen atoms, Vpx,py is also small 



(0.1%), indicating the orbitals are equally populated (data 
not shown). Vpx,pz and Vpy^pz have smaller values (1%) 
compared to the apical oxygen atoms in all distortion 
ranges examined. For large, physically unreasonable dis- 
tortions, the Pz orbital not participating in the in-plane 
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FIG. 4. Orbital polarizations for the apical oxygen atoms Vpx.pz (a) and (b), Vpy.pz (c) and (d), and Fe(l) site 'Pd2.do (e) 
and (f) for breathing distortions (BD) coupled with in-phase (upper panels) and out-of-phase (lower panels) rotation patterns, 
respectively. aPaPc' rotations induce charge transfer to the p^ orbital aligned along the Fe-0 apical bond axis; Fe eg polarization 
is weakly affected by rotation sense. Filled (red) circles indicate the BD amplitude and the rotation angles for the structures 
whose DOS are given in Fig. 2(b) and (c). 



Fe-O bond eventually becomes populated. Note that 
oxygen orbital polarizations involving pz are not affected 
by small breathing distortion values; however, coupling 
a large breathing amplitude with large rotation angles 
makes V maximum. 

We now examine the Fe 3d-orbitals with Cg symme- 
try. Fe(l) orbital polarization 7^d2,do[Fe(l)] increases with 
breathing distortion [Fig. 4(e-f)], i.e., charges flow from 
the apical to equatorial bond axes. With decreasing ampli- 
tude of the rotation angle Oxy^ the d^^'^-r^ becomes pref- 
erentially occupied [Fig. 3(d)]. For Fe(2) the 75^2,^0 [Fe(2)] 
decreases as the charge flow reverses from the equatorial 
to apical bond. Both polarizations are minimal for small 
rotation angle values (further from cubic symmetry), yet 
the differential occupation appears to be weakly depen- 
dent on the rotation sense. The out-of-phase rotations, 
however, influence the eg orbitals more than the in-phase 
ones. 

To summarize to here, our analysis reveals that the 
apical oxygen atoms lying along the rotation axis are 
mainly sensitive to the octahedral rotation sense, with 
out-of-phase rotations favoring charge flow to the Fe-0 
apical bond. The orbital occupancy of the equatorial 
oxygen atoms is weakly affected by the rotation pattern 
sense. '" Moreover, the iron sites are largely sensitive to 
the breathing distortion, whereby the more expanded an 
octahedron becomes, the greater is the charge transfer 
to the apical bond axes. This effect increases the orbital 
polarization of the corresponding Fe cations; it can also 
be more readily achieved with the a^a^c~ rather than 
the a^a^c"*" rotation pattern. These results suggest that 
sense of the octahedral rotation mainly affects the oxygen 



atoms lying along the rotation axis that "screen" the 
metal cations, which are unable to "see" the rotation 
sense effects. 



B. Orbital Polarization with Multiple Rotations 



In the experimental CaFeOa structure, the BD is 
present together with the a~ a~ c'^ rotation pattern, where 
both rotation senses are present and out-of-phase rota- 
tions are found about two Cartesian directions. We now 
analyze the orbital polarization in CFO structures where 
the BD is coupled with the a^a^c^^ rotation to evalu- 
ate the conclusions about the different rotation sense 
discussed above in the presence of multiple rotation axes. 

We find polarizations Vpx.py, Vpx.pz and Vpy^pz of the 
apical oxygen atoms increase when the rotation angle 9z 
approaches 180° (data not shown). '^ As before, the out- 
of-phase rotations will transfer charge from the orbitals 
lying in the plane normal to the rotation axes to the 
orbital lying along the axes. In this case, the rotations 
in the xy-plane induce a flow of charge from pz to both 
the px and Py orbitals, i.e. from the Fe-0 apical bond 
into the plane orthogonal to it and containing the apical 
oxygen atom and its Px and Py orbitals. On the other 
hand, the equatorial oxygen atoms experience charge flow 
to their corresponding Fe-0 bonds. Concerning the two 
Fe sites, decreasing the 6z rotation angle moves charge 
from the apical (d^z^-r^ orbital) to the equatorial ((i^2_j,2) 
axis. Here, the BD has the role of enhancing the charge 
transfer induced by the a~a~c'^ rotation pattern. This 
last result, together with the analysis of the single in- 



phase and out-of-phase case discussed above, reveals that 
the out-of-phase rotation pattern favors a flow of charge 
towards its corresponding rotation axis. The rotation 
"sense" becomes an important structural feature over the 
orbital population. 

In CaFeOa, the Fe-0 equatorial bonds are known 
to play a fundamental role in the metal-insulator 
transition:^'' a change in the equatorial bond-lengths 
together with the octahedral rotations affects the 
dx^ -y2 / (px ,Py) orbital hybridization responsible for the 
formation of a highly covalent electronic configuration, 
gapping the system. An alternative route for tuning the 
band gap emerges from our results. If the Fe-0 bond 
lengths are not altered, then the a''a~c^ out-of-phase 
rotation pattern controls the amount of charge filling 
the equatorial orbitals, hence the degree of hybridization, 
which determines the stability of the insulating phase. 



C. Strain— rotation and Orbital Polarization 
Coupling 

Strain can enhance the octahedral rotations in 
perovskites,'" which in turn alter the band structure.'^'"''*' 
Tensile bi-axial strain promotes the a~a~c'^ rotation pat- 
tern in c-oriented orthorhombic perovskites thin films. " 
To explore the effect of strain on rotation-induced V, we 
compute the equilibrium structure of CaFeOs under biax- 
ial tensile strain, relaxing ionic positions and the c-axis 
length. 

First, we evaluate the effect of strain on the atomic 
structure. We find that 6z deviates more from 180° with 
increasing tensile strain [Fig. 5(a)]. Consistent with our 
previous observation, when the BD and a^a~c^ rotations 
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FIG. 5. BD amplitudes and 9^ angles (a), orbital polarization 
(%) of (b) Fe(l) site, (c) equatorial (eq.) and (d) apical (ap.) 
oxygen atoms for relaxed CaFeOs structures under biaxial 
strain. The shaded area indicates strain values for insulating 
structures; the MIT occurs between 1-1.5% strain. 



are active, the orbital polarization is enhanced along 
the Fe-0 equatorial bonds with greater tensile strain. 
Importantly, we note that within the insulating phase as 
strain is increased, the BD largely remains constant while 
the octahedral rotations are enhanced to accommodate 
the strain energy [Fig. 5(a)]. The result of the decreased 9z 
rotation angle is further charge localization in the d^i_yi 
[Fig. 5(b)] and Px and Py orbitals of equatorial oxygen 
atoms [Fig. 5(c)], hence about the equatorial Fe-0 bonds 
supporting the insulating state and leaving unaltered the 
polarization of apical oxygen atoms [Fig. 5(d)]. 

Our results show that the O 2p and Fe Cg orbital po- 
larizations are highly sensitive to the strain state and 
are mediated largely by the change in the out-of-phase 
octahedral rotations {9z)- The tensile strain, enhances 
the a~a~c^ rotation pattern, induces a charge transfer 
into those orbitals responsible for the gap (i.e. dx2-y2, 
and Px/Py of equatorial oxygen atoms): beyond a critical 
value, the charge becomes delocalized and the insulating 
gap closes to form a metallic state. We observe that the 
insulating phase is lost for strain values greater than 1%, 
which is in contrast to that seen in rare-earth nickelates,^'' 
suggesting the rotation-orbital polarization interaction 
and the BD may be either complementary or antagonistic 
to the MIT. 

The out-of-phase rotations are also responsible for the 
broadening of the Px and Py orbitals about Ep , which we 
conjecture would reduce the sharpness of the first-order 
phase transition. Note that the critical transition temper- 
ature for the MIT is also believed to be correlated with the 
size of the BD, our results therefore suggest that strain- 
engineering of the MIT temperature could be challenging 
in ferrates. On the other hand, the structural transition 
temperature can be kept fixed with increasing strain while 
the electronic transition can still be tuned by means of 
the octahedral rotation-induced orbital polarizations. 



IV. CONCLUSIONS 

In summary, we have shown how the sense of the oc- 
tahedral rotation pattern and not only the magnitude of 
the rotation distortion is important in engineering the low 
energy electronic band structure. We find that while the 
in-phase rotations weakly affect the local electronic dis- 
tribution, the out-of-phase rotation pattern directs a flow 
of charge towards its corresponding rotation axis. The 
size of the rotation-sense orbital polarization we compute 
should be experimentally detectable with orbital reflec- 
tometry techniques,''* across the orthorhombic {a^a^C^) 
to monoclinic {a~b~c'^) phase transition of CaFeOs near 
room temperature. Here, the out-of-phase rotation pat- 
tern is no longer uniform along [110]. For charge transfer 
perovskite oxides, control over the phase of the rotations 
via heteroepitaxial strain, substrate proximity effects, or 
superlattice formation provides a promising alternative 
route to tailor the local orbital polarizations, hence elec- 
tronic transitions, in correlated materials near electronic 



phase boundaries. Due to the complexity of these inter- 
actions driven by cation substitution, we hope this work 
motivates further experimental studies to evaluate the 
mechanisms proposed here, and the degree to which they 
compete or cooperate across a MIT. 
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